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LONG-TERM  GOALS 

The  long  term  goal  of  this  proposed  research  are  twofold:  (1)  develop  an  advanced  multi-physies 
model  with  a  multi-numerieal  solution  teehniques  approach  to  predict  nonlinear  dynamie  behavior 
of  impaet  burial  and  flow-indueed  motion  of  flexible  struetures  (mines)  and  surrounding  sediments 
(sand)  in  the  marine  environment;  and  (2)  ealibrate  resulting  models  with  experimental  and  field 
measurements.  The  predietive  eapability  developed  in  this  researeh  will  eventually  be  integrated 
into  an  overarehing  eomputational  framework  for  the  analysis  and  simulation  of  the  dynamie 
behavior  of  naval  systems  in  the  marine  environment  of  arbitrary  water  depth. 

OBJECTIVES 

The  ability  to  detect  and  neutralize  bottom  mines  is  critieally  dependent  on  the  scour  pattern  around 
the  mine  and  on  the  degree  of  the  mine  burial.  The  objeetives  of  the  first  two-years  of  this  researeh 
are  to  first  gain  an  understanding  of  the  strueture-fluid-sediment-seabed  interaction  phenomenon  by 
condueting  a  literature  survey  to  determine  state-of-the-art  analytieal  and  numerical  modeling 
eapabilities.  Then  evaluate  the  eurrent  status  of  analysis  and  simulation  software  LS-DYNA  in 
terms  of  its  predietion  eapabilities  for  eoupled  dynamie  motions  of  submerged  mines  on  a  seabed  in 
the  marine  environment  and  to  identify  immediately  needed  developments  for  further  improvement. 
Simultaneously,  survey  and  examine  laboratory  and  field  experiments  of  the  nonlinear  dynamic 
behavior  of  impact  burial  and  flow-induced  motion  of  mines  and  scour  of  surrounding  sediments  in 
the  marine  environment. 

APPROACH 

The  dynamie  behavior  of  submerged  mines  and  their  surrounding  sediments  has  been  of  interest  to 
the  Navy  in  recent  years  (Inman  and  Jenkins,  2002).  During  1960s-1980s,  U.S.  Naval  seientists 
developed  three  burial  predietion  models  that  can  be  used  for  mine  eountermeasure  tactical  planning 
and  for  development  of  environmental  support:  impact  burial,  sand  ridge  migration,  and  wave- 
induced  scour  (Richardson  and  Briggs  2000).  Capabilities  for  aeeurate  modeling  and  prediction  of 
motions  of  mines  and  their  effects  on  the  overall  burial  behavior  at  the  sea  bottom  eonsidering  the 
detailed  physics  of  the  entire  coupled  mine-fiuid-sediment-seabed  “system”  is  needed.  The  eode 
LS-DYNA,  which  is  an  advanced  nonlinear  finite-element  based  eommereially  available  numerieal 
software,  has  arguably  the  best  solid  meehanies,  eontact  and  impaet  models  among  eommereially 
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available  software,  and  has  growing  fluid  mechanie  features.  It  also  has  an  exeellent  eomputational 
framework,  and  pre-  and  post-proeessing  eapabilities  for  us  to  take  advantage  of  in  developing 
additional  analysis  eapabilities  to  address  our  researeh  and  applieation  needs. 

The  overall  objeetives  of  this  researeh  are  to  be  aehieved  through  the  following  tasks:  (1)  Conduet  a 
thorough  review  of  literature  on  mine-fluid-sediment-seabed  interaetion  to  determine  the  state-of- 
the-art  analytieal  and  numerieal  modeling  eapabilities;  (2)  Evaluate  the  eurrent  predictive  capability 
of  LS-DYNA  for  coupled  dynamic  motions  of  submerged  mines  on  a  seabed  in  the  marine 
environment;  (3)  Identify  development  needs  of  LS-DYNA  for  coupled  fluid-structure-seabed 
interaction  (mine  burial  and  motion  prediction)  applications;  (4)  Plan  and  implement  the  immediate 
needed  developments  within  a  two-year  timeframe;  (5)  Examine  existing  laboratory  and  field 
experiments  available  from  ONR  on  motions  of  mines  on  seabed  for  numerical  model  calibration; 
(6)  Compare  numerical  predictions  of  resulting  numerical  modules  developed  with  laboratory 
experiment  and  field  data  to  validate  and  calibrate  the  numerical  models  for  further  evaluation;  (7) 
Conduct  parametric  study  of  characteristic  burial  and  exposure  times,  and  other  parameters  as 
appropriate;  and  (8)  Analyze  and  document  research  results  and  develop  a  future  research  plan  for 
further  needed  modeling  improvements. 

WORK  COMPLETED 

In  the  first  two  years  of  this  project,  we  have  focused  on  evaluation  of  the  current  status  and 
development  of  immediately  needed  capabilities  of  a  multi-physics  code  LS-DYNA  for  impact 
burial  and  flow-induced  motion  predictions  of  mines  and  its  surrounding  sediments.  While  the  study 
of  impact  burial  has  been  quite  active  and  successful  from  both  analytical  and  experimental 
perspectives,  the  study  of  motions  and  scour  of  buried  mines  has  been  restricted  mostly  to  empirical 
and  experimental  perspectives.  We  conducted  a  review  of  the  literature  on  impact  burial  and  mine- 
fluid-sediment-seabed  interaction  paying  special  attention  to  software  developed  by  researches  as  it 
is  needed  to  analyze  the  coupled-ESI  interaction  problem.  We  also  examined  the  LS-DYNA 
software  to  evaluate  its  current  capabilities  from  the  prospective  of  modeling  the  dynamics  of  the 
coupled  LSI  system.  There  are  several  soil  material  models  in  LS-DYNA,  which  are  currently  under 
investigation,  to  determine  if  one  of  these,  if  enhanced,  would  meet  the  requirement  of  a 
cohesionless  soil.  A  brief  list  of  all  the  soil  models  in  LS-DYNA  is  given  in  the  appendix.  The  three 
common  laboratory  tests  that  are  performed  to  get  the  various  material  model  parameters  for  the 
constitutive  equations  are  (i)  hydrostatic  compression  (ii)  triaxial  compression/extension  and  (iii) 
uniaxial  strain.  The  material  model  parameters  that  can  be  calibrated  to  these  data  are  used  in  the 
LS-DYNA  soil  constitutive  models.  A  selected  number  of  these  models  were  reviewed  as  possible 
soil  model  candidates,  from  the  simplest  (material  model  5)  to  the  most  complex  (material  model 
25).  Most  of  the  candidate  models  in  LS-DYNA  are  extensions  of  material  model  5  (soil  and  foam). 
Two  of  the  exceptions  to  this  are  model  16  (pseudo-tensor  model)  and  model  25  (geological  cap 
model).  Model  5  and  its  extensions,  model  14  (soil  and  foam  with  failure),  and  model  79  (hysteretic 
soil)  are  basically  an  analytical  pressure  dependent  yield  surface.  Both  of  these  models  were  indeed 
found  to  be  unstable  in  unconfined  states.  Material  model  14  and  model  78  (soil  and  concrete)  are 
for  analyses  in  high-pressure  regions.  Model  14  was  evaluated  previously  and  it  was  found  that  it 
does  not  simulate  low/zero-pressure  behavior  well.  Material  model  25  (geological  cap  model)  is  a 
complex  model,  but  does  handle  low-confinement  behavior  well.  Model  25  cannot  simulate  strain 
softening.  Despite  this  shortcoming,  it  is  the  only  existing  model  that  could  simulate  much  of  the 
basic  behavior  needed  for  the  present  purpose.  Because  of  the  cap/shear  surface  intersection  (non¬ 
smooth)  in  model  25,  there  are  many  corners  in  the  yield  surface.  These  comers  cause  the  algorithm 
to  be  very  complex  and  inefficient.  Although  there  are  several  other  soil  material  models  existing  in 
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LS-DYNA,  only  those  whose  material  characterization  is  readily  available  were  investigated.  Based 
on  these  observations,  it  was  concluded  that  none  of  the  existing  models  are  ideal  for  modeling 
cohesionless  sand,  and  it  was  decided  to  test  a  few  of  these  models  for  the  mine  burial  tests  to 
determine  the  capability  and  needed  improvements  of  LS-DYNA  to  model  cohesionless  sand.  In 
order  to  employ  these  models  for  mine-sand  interaction  simulations,  experimental  results  to 
determine  the  material  characterization  are  needed.  We  assisted  in  an  experiment  on  scour  of  sandy 
seabed  around  a  vertical  cylinder  and  have  access  to  data.  We  are  also  in  communication  with  the  PI 
of  an  on-going  field  experiment  on  mine  scour  being  conducted.  As  a  pure  continuum  approach  is 
not  able  to  capture  the  complex  motion  of  the  sand  particles  during  the  settling  and  burial  of  mines, 
we  are  currently  focusing  our  study  on  a  discrete  approach  based  on  the  smoothed  particle 
hydrodynamics  method  for  the  mine  burial  problem. 

RESULTS 

Understanding  the  burial  and  re-exposure  of  submerged/bottom  mines  placed  in  nearshore  waters,  is 
of  paramount  interest  to  the  navy.  The  ability  to  detect  and  neutralize  bottom  mines  is  critically 
dependent  on  the  scour  pattern  around  the  mine  and  on  the  degree  of  the  mine  burial.  In  order  to 
understand  the  extent  to  which  LS-DYNA  can  model  scour  scenario  and  at  the  same  time 
comprehend  the  capability  of  the  software  to  handle  various  soil  models,  a  wave  generation 
experiment  was  conducted  to  simulate  a  simple  experiment.  The  model  contains  the  mine,  fluid  and 
a  piston  type  wave  generator  which  is  defined  in  rigid  shell  elements.  Geometry  and  dimension  of 
the  wave  generating  tank  was  so  chosen  to  generate  a  1 .20m  high  wave  from  the  sea  level  (Sunao 
Tokura,  2005).  The  penalty  coupling  was  used  to  define  the  fluid-structure  interaction.  The  water 
was  modeled  by  using  solid  brick  elements.  A  void  domain  was  created  on  top  of  the  water  domain 
to  capture  the  wave  elevation.  The  water-void  mesh  was  modeled  as  a  rectangular  mesh.  The  length 
of  the  wave  generating  tank  was  40.0m  and  a  width  of  6.0m.  The  height  of  the  water  column  was 
2.0m  and  size  of  the  void  domain  was  3.0m.  The  mine  was  rectangular  in  shape  (4.0m  x  0.8m).  The 
depth  of  the  soil  was  1  .Om  and  the  mine  was  allowed  to  rest  on  the  soil  and  gravity  loading  was 
applied  to  allow  for  the  mine  to  settle,  before  running  the  wave  maker.  The  numerical  wave  tank  has 
a  non-refiective  boundary  condition  at  the  far  end  of  the  solution  domain  in  order  to  absorb  the 
wave  energy  and  prevent  reflection.  The  other  boundaries  of  the  solution  domain  are  along  the  side 
walls  of  the  tank  (impermeable  no-slip  boundary).  Figure  1  shows  the  computational  mesh.  The 
tests  were  repeated  for  three  soil  models,  however,  as  the  other  two  tests  produced  results  that  are 
similar  in  nature,  only  the  results  from  the  first  soil  model  is  shown.  Fluid  density  animation  plots 
for  the  Soil  and  Foam  model  at  various  time  steps  are  shown  in  Figure  2.  The  existing  models  are 
based  on  the  Drucker-Prager  soil  models  which  are  directly  based  on  either  the  Mohr-Coulomb 
failure  surface  or  the  Modified-Mohr-Coulomb  failure  surface  which  cannot  be  used  to  properly 
understand  the  saturated  sand  behavior.  As  it  is  difficult  to  extrapolate  dry  sand  material  properties 
to  those  for  fully  saturated  sand,  the  option  of  treating  saturated  sand  with  the  properties  of  water 
(low  shear  strength  and  high  density)  was  also  tested.  The  simulations  thus  far  utilizes  a  continuum 
approach  to  model  each  of  the  fluid,  sand,  solid  mediums  using  an  Arbitrary  Lagrangian  and 
Eulerian  (ALE)  approach.  As  the  use  of  ALE  techniques  is  not  sufficient  to  fully  analyze  the  scour 
around  the  mine,  a  robust  soil  model  is  needed  to  fully  capture  the  scour  scenario.  Apart  from  this 
an  additional  material  characterization  is  needed  to  fully  understand  the  capability  of  ES-DYNA  in 
handling  saturated  sand  behavior.  The  key  for  these  kinds  of  modeling  studies  is  a  soil  model  that 
includes  the  pore  pressure  effects  in  the  saturated  sand  and  clearly  none  of  the  existing  ES-DYNA 
geomaterial  models  include  pore  pressure  into  its  effect.  However,  it  has  been  found  that  a  pure 
continuum  approach  is  not  sufficient  to  capture  the  complex  motion  of  the  sand  particles  during  the 
settling  and  burial  of  mines.  A  discrete  approach  such  as  smooth  particle  hydrodynamics  (SPH), 
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which  relies  on  the  equations  of  state  to  represent  the  material  behavior,  is  being  investigated  as  a 
viable  alternative.  Retaining  a  similar  wave  generation  experiment  over  a  sand  bed,  a  eomparative 
study  to  understand  seour  seenario  around  a  solid  objeet  (mine  burial)  will  be  eondueted  using  the 
SPH  method  available  in  LS-DYNA. 

SPH  is  a  purely  Lagrangian  method  developed  during  the  1970’s  (Lucy  1977;  Gingold  and 
Monaghan  1977)  in  astrophysics  to  study  the  eollision  of  galaxies  and  the  impaet  of  bolides  on 
planets.  Unlike  finite  element  methods,  the  SPH  method  discretizes  the  deformable  objeet  with 
partieles  rather  than  connected  meshes.  The  main  advantage,  however,  arises  direetly  from  its 
Lagrangian  nature,  sinee  sueh  an  approaeh  ean  tackle  diffieulties  related  with  laek  of  symmetry, 
large  voids  that  may  develop  in  the  field,  and  a  free  water  surfaee  mueh  more  effieiently  than 
Eulerian  methods  ean.  The  SPH  teehnique  available  in  LS-DYNA  has  been  applied  for 
ineompressible  dynamic  fluid  flow  patterns.  A  series  of  numerieal  tests  have  been  carried  out  to 
examine  the  ability  and  effieieney  of  SPH  formulations  in  simulating  fluid  dynamie  problems.  The 
SPH  method  with  various  formulations  ean  simulate  different  dynamie  fluid  flow  problems,  sueh  as 
inviseid  or  viseous  flows,  compressible  or  incompressible  flows.  The  SPH  option  in  LS-DYNA  was 
used  to  simulate  the  Poiseuille  flow  and  Couette  flow.  The  SPH  option  in  LS-DYNA  demands  an 
equation  of  state  and  the  termination  time  was  set  to  0.5  seeonds.  In  the  elassieal  hydrodynamies, 
the  flow  veloeity  at  a  point  in  the  Poiseuille  flow  or  Couette  ean  be  obtained  by  solving  the  Navier- 
Stokes  momentum  equation.  In  this  applieation  the  cubie  spline  funetion  is  used.  The  Poiseuille 
flow  involves  flow  between  two  parallel  stationary  infinite  plates  placed  at  y  =  0  andy  =  / .  The 
initially  stationary  fluid  is  driven  by  a  body  foree  F  (e.g.,  pressure  differenee  or  external  foree), 
gradually  flows  between  the  two  plates,  and  finally  arrives  at  a  steady  state.  Couette  flow  is  another 
often-used  benehmark  CFD  problem.  It  is  a  flow  between  two  initially  stationary  plates  plaeed  at 
y  =  0  and  y  =  /  when  the  upper  plate  moves  at  a  eertain  eonstant  veloeity  (  Vq).  The  simulations  of 
Poiseuille  and  Couette  flow  show  that  this  approaeh  ean  be  furthered  to  understand  the  seour  around 
a  bottom  mine. 

IMPACT/APPLICATIONS 

The  advaneed,  state-of-the-art  FE  eode  FS-DYNA  adopted  in  this  projeet,  when  fully  developed, 
will  enhanee  the  modeling,  predietion,  operation  and  eontrol  eapabilities  of  the  eomplex  mine-fluid- 
sediment-seabed  interaetion  in  general  and  the  numerieal  simulations  of  fiow-indueed  mine  burial 
impaets  in  partieular.  The  3-D  numerical  codes  being  developed  will  provide  additional  tools  to 
ealibrate  and  validate  the  accuracy  of  the  numerieal  predietions  of  the  modules  with  laboratory 
experiment  and  field  data.  Figure  3  shows  the  eomputational  model  for  eoupled  fiuid-strueture 
interaetion  (FSI)  problem  with  respeet  to  the  mine  burial  in  sand  sediments.  From  the  eomputational 
sketch  it  can  be  seen  that  the  fluid  domain  followed  by  the  strueture  and  the  sand  ean  be  modeled 
using  four  different  eomputational  methods.  The  fluid  domain  will  be  modeled  using  the  fully 
nonlinear  potential  flow  (FNPF)  or  the  boundary-element  method  (BEM).  Reynolds  averaged 
Navier-Stokes  (RANS)  and  the  partiele  finite  element  method  (PFEM)  will  be  used  in  the 
water/mine/sand  domain.  Sand  and  the  geomaterials  around  the  sand  will  be  modeled  using  the 
Smoothed  Partiele  Hydrodynamics  method  (SPH).  The  resulting  numerieal  predictive  eapability 
will  be  ineorporated  into  an  overarehing  eomputational  framework  for  the  analysis  and  simulation 
of  the  dynamie  behavior  of  naval  systems  in  the  marine  environment  of  arbitrary  water  depth. 
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TRANSITIONS 


Analysis  and  simulation  capabilities  developed  in  this  research  can  be  useful  to  the  various  units  of 
the  Navy  pertinent  to  mine  deployment,  detection,  burial  clearance  proeess  studies. 

RELATED  PROJECTS 

The  analysis  and  simulation  eapabilities  developed  in  this  researeh  will  be  ineorporated  into  a 
eompanion  project  (N00014-07-1-0008)  on  the  development  of  an  overarehing  eomputational 
framework  for  analysis  and  predietion  of  dynamic  motions  of  naval  systems  in  the  marine 
environment  in  arbitrary  water  depth. 
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Fig.  1  The  computation  domain  (Isometric  View) 
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Fig.  2  Fluid  density  animation  plots  for  wave  propagation  over 
a  bottom  mine  at  various  time  steps  (Soil  and  Foam  model) 
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PANS:  Reynolds  Averaged  Navier  Stokes 


FEM:  Finite  Element  Method 


Fig.  3  Computation  model  (various  methods  that  will  be  used)  for  the  coupled  Fluid  Structure  Interaction  (FSI)  problem 
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